DIRAC AND MAGNETIC SCHRODINGER OPERATORS ON 

FRACTALS 



MICHAEL HINZi AND ALEXANDER TEPLYAEV^ 



Abstract. In this paper we define (local) Dirac operators and magnetic Schrodinger 
Hamiltonians on fractals and prove their (essential) self-adjointness. To do so we use 
the concept of 1-forms and derivations associated with Dirichlet forms as introduced 
by Cipriani and Sauvageot, and further studied by the authors jointly with Rockner, 
lonescu and Rogers. For simplicity our definitions and results are formulated for the 
Sierpinski gasket with its standard self-similar energy form. We point out how they may 
be generalized to other spaces, such as the classical Sierpinski carpet. 
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1. Introduction 

The aim of the present paper is to introduce natural (local) Dirac and magnetic 
Schrodinger operators on fractal spaces and to prove that they are (essentially) self- 
adjoint. Our analysis uses the concept of 1-forms in the content of the Dirichlet form 
theory, and is based on recent results on 1-forms and vector fields |T8l HH [201 EO] and |17] , 
respectively. 

To make the paper more accessible and to approach the most interesting classical exam- 
ples, we formulate our definitions and results for the Sierpinski gasket, and later provide 
some remarks how to modify them for more general fractals and other spaces. This is par- 
ticularly straightforward to do for species with a resistance form, in the sense of Kigami 
[MIESIISS] (see also SSI EB] ) , such as the classical two-dimensional Sierpinski carpet, 
but many results are valid for much more general spaces. In particular, extending our 
results for spaces that are not locally compact will be subject of future work. 

Our space of 1-forms is a Hilbert space, which allows to identify 1-forms and vector 
fields, and to introduce other notions of vector analysis, as recently done in |17] (which 
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generalizes earlier approaches to vector analysis on fractals, see [521 [551 EHl [921 [98]). This 
is a part of a comprehensive program to study vector equations on general non-smooth 
spaces which carry a diffusion process or, equivalently, a local regular Dirichlet form. 

The study of the Laplacian on fractal graphs was originated in physics literature (see 
[51 [TH [2S1 [SSI ES])) cind for a selection of recent mathematical physics results see 
[21 [21 [221 [23 [ini [511 and references therein]. Among the problems where fractal spaces 
seem to appear naturally we would like to mention, in particular, the spaces of fractional 
dimension appearing in quantum gravity [6l [691 [821 and references therein]. Besides that, 
our motivation is coming from the theory of quantum graphs [301 [211 1251 [121 [591 EOl [SH [621 
[63l [Ml [761 and references therein]; from the spectral theory on fractals [H [23 [HI [571 Ell 
[671 [721 [731 [751 EH [851 and references therein]; form some questions of non-commutative 
analysis [221 [23 [HI [13 [201 [501 and references therein] and the theory of spectral zeta 
functions [23 [SSI [SOI [SZ]; and from the localization problems [H [73 113 ESI [13 and 
references therein]. 

Recall that roughly speaking, the Dirac operator is defined as the square root of the 
Laplace operator. (Note, however, that classical Dirac operator for diffusions is a local 
operator, which excludes the possibility of using the spectral theorem to define it.) De- 
pending on context and purpose it appears in various formulations with possibly different 
complexity and sign conventions. On the real line D = —id/dx may for instance be 
regarded as the Dirac operator. Given a Riemannian manifold M, its tangent bundle 
AT*M can be turned into a Clifford module, and the associated Dirac operator is defined 
as D = d + d*, where d is the exterior derivative and d* is its adjoint, cf. [11]. For a 
spin 1/2-particle in the plane the Dirac operator is given hy D = —ia^d/dxi —iayd/dx2, 
where and ay are the respective Pauli spin matrices, [21]. More generally, it may be 
defined for spinor bundles over spin manifolds, see [13 [21] or [28] for background and 
details. 



Dirac operators on discrete graphs have for instance been considered in p5], Section 

4] with a strong emphasis on connections to noncommutative geometry. The paper [8Tj 
follows a similar spirit and considers related spectral triples and Connes metrics. More 
recently Dirac operators on discrete graphs and related index theorems have been studied 
in [7^. In this reference they act on a tensor product of form Hq® Hi, where Hq and Hi 
are Hilbert spaces of functions on the vertices and edges, respectively. Roughly speaking, 
the discrete difference operator d : Hq ^ Hi plays the role of the exterior derivative. 
Denoting its adjoint by d* : Hi ^ Hq , the associated Dirac operator is then defined on 
Ho © Hi by 



and as a consequence yields the matrix Laplacian acting on Hq® Hi. In a somewhat 
similar fashion [76] also investigates Dirac operators and index theorems on quantum 
graphs (often referred to as metric graphs or quantum wires, [53 [63 [621 [S2]), now within 
the context of suitable Sobolev spaces. A preceding article dealing with index theorems 
on quantum graphs is , and a related much earlier reference for Dirac operators is [T7] . 
Different quantization schemes are reviewed in [12]. 
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There is an extensive literature on magnetic Schrodinger operators. In the simplest 
cases (such as for bounded and sufficiently integrable potentials A and V) the essential 
self-adjointness of magnetic Schrodinger Hamiltonians 



on Euclidean space can be deduced from classical perturbation theorems in Hilbert spaces, 
cf. [Sol Section X.2]. More sophisticated pointwise methods can be found in |80( Section 
X.3 and X.4]. Essential self-adjointness results for operators on manifolds may be found 
in the comprehensive paper [89], see also [39] for singular potentials. Discrete mag- 
netic Schrodinger operators on lattices and graphs have for instance been discussed in 
[HI ttH [221 12S1 SSI EHl MM ■ The paper [SS] introduces discrete magnetic Laplacians on the 
two-dimensional integer lattice, proves that they have no point spectrum and compares 
them to almost Mathieu operators and to one-dimensional quasi-periodic operators. This 
is closely connected to the (one-dimensional) ten martini problem, solved in Peri- 
odic magnetic Schrodinger operators on the two-dimensional integer lattice are treated 
in [91] and their spectra, typically of band or Cantor type, are studied. In [13] magnetic 
Schrodinger operators on graphs are considered. Under some conditions the analyticity 
of the bottom of their spectra is verified and relations to corresponding operators on a 
quotient graph (by a suitable automorphism group) are discussed. The paper [2n] also 
investigates discrete magnetic Laplacians and Schrodinger operators on graphs, compares 
their spectra and heat kernels to the original graph Laplacians, defines related Novikov- 
Shubin invariants and establishes a long term decay result for the heat kernel trace. In 
[22] the essential self-adjointness of a discrete version of (|2|) is shown, based on a previous 
result [21] for operators with zero magnetic potential. Reference [22] also discusses gauge 
invariance in terms of holonomy maps. First steps towards magnetic Schrodinger oper- 
ators on fractals had been taken in [9j and ^Jj by studying them on infinite Sierpinski 
lattices. Some decimation techniques for the spectrum and related numerical experiments 
can be found in [9] . The paper [TT] sets up a renormalization group equation for the mag- 
netic Laplacian and discusses relations to superconductivity. Another branch of literature 
concerns quantum graphs, see [SHI EHl E21 ES] • The paper [M] introduces magnetic Lapla- 
cians on metric graphs and, based on results in [60], provides a matrix criterion for the 
boundary conditions to characterize self-adjointness. In [16] a metric graph point of view 
is used to provide a comprehensive study of the two-dimensional periodic square graph 
lattice with magnetic fields. The paper [30] shows that any self-adjoint vertex coupling 
on a metric graph can be approximated by a sequence of magnetic Schrodinger operators 
on a network of shrinking tubular neighborhoods. 

For prototype examples of fractal sets carrying a diffusion not even the forms of a Dirac 
operator and a magnetic Laplacian had been clear. This is due to the fact that Laplace 
operators had been studied in several papers and books (see for instance [TUl 1521 ESI [SH 
mi [91] and the references therein), but definitions and results concerning analogs of first 
order differential operators (gradients) were sparse [521 [SSl [921 [M], and hardly flexible 
enough to fit a sufficient functional analytic context. 

In [19] and [2Q] differential 1-forms and derivations based on Dirichlet forms had been 
introduced. In these papers a Hilbert space "H of 1-forms is constructed as, roughly speak- 
ing, the completion of the tensor product oi the space J-" of energy finite functions, 
a concept that leads to an L2-theory, see for instance [IS] for further explanations. This 
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approach has been studied further in [HI [50] and also in |17] , where related notions of 
vector analysis are proposed. In this context the desired objects can be defined. More 
precisely, our Theorems 13.11 and 14.11 below state that under some conditions, analogs of 
([T]) and ([2]) define essentially self-adjoint operators on suitable Hilbert spaces of functions 
and vector fields on fractals, respectively. 

In the next section we review the approach of [T^ |2D] to 1-forms based on energy for 
the specific example of the Sierpinski gasket K with its standard self-similar energy form 
{S,J^). We recall the definitions of the gradient and divergence operators from |17j and 
the energy Laplacian for functions. In several places we provide auxiliary formulations in 
terms of harmonic coordinates. In Section [3] we define a related Dirac operator D that acts 
on the tensor product L2,£{K, v) ® Tic of the spaces of complex- valued square integrable 
functions (with respect to the Kusuoka measure z/) and complex-valued vector fields on 
K. Theorem 13.11 proves it it is self-adjoint. In Section H] we first provide a priori esti- 
mates necessary to introduce a bilinear form S"-'^ associated with a magnetic Schrodinger 
Hamiltonian H"^'^ on K. Then we establish a result on its essential self-adjointness on 
L2,c{K, I'), Theorem 14. H which merely follows from our definitions, preliminary estimates 
and a simple KLMN theorem. Finally we prove sort of a gauge invariance result. Theo- 
rem 14. 2[ Section [5] contains some instructions how to generalize the presented results to 
arbitrary finitely ramified fractals carrying a regular resistance form. 

In this paper we generally intend to provide a basic setup to study Dirac operators and 
magnetic fields on fractals. We do not discuss questions regarding the spectrum, refined 
pointwise statements or approximations. These topics will be addressed elsewhere. 

To simplify notation, sequences or families indexed by the naturals will be written 
with index set suppressed, e.g. (a„)„ stands for (a„)„gN- Similarly, lim„ a„ abbreviates 

2. Vector analysis on the Sierpinski gasket 

This section recalls a few items of the concept of 1-forms and vector fields based on 
Dirichlet forms as studied in [181 ISl 1201 EO] and [17], respectively. For simplicity we 
formulate definitions and results for the Sierpinski gasket, some comments on finitely 
ramified fractals are provided in Section O For investigations of other physical models on 
the Sierpinski gasket see for instance [321 [23] • 

Let {pi,P2,P3} be the vertex set of an equilateral triangle in M^. The Sierpinski gasket 
K is the unique nonempty compact subset of such satisfying the self-similarity relation 

3 
i=l 

where ipi{x) = x/2 + pi/2. For our purposes the embedding in ]R2 is 

inessential, only 

the associated post critically finite self-similar structure {K, {1, 2, 3} , f2, ^3}) in the 
sense of [53j matters. Let {S,J-') denote the standard self-similar energy form on K, 
obtained as the increasing limit 



£{u) := \im £n{u) 

n 
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of a sequence of rescaled discrete energies 

on approximating graphs. For precise definitions and background we refer to [13l |20l [521 
[55] . The form {£,J-') is a regular resistance form in the sense of [SS]. In particular, E^^"^ 
is a Hilbert norm on the space J-"/ ~ obtained from J-' by factoring out constants, and 
there is some constant c > such that 

(3) ll/lloo<c^(/)'/' 

for all / G J-"/ ~. The space J-" is a dense subalgebra of C{K), and in particular 

(4) Sifgy/' < £{ff'^ hh^^K,.) + \\f\\L^iK,.)m"\ f,9^J'- 

For any / G J-" we can define a unique (nonatomic) Borel energy measure vj on K, see 
for instance [96] , and polarization yields mutual energy measures z//^g for f,g^J^. 

The space of nonconstant harmonic functions on K is two dimensional. Let {hi, /12} be 
a complete energy orthonormal system for it. The Kusuoka energy measure v is defined 
by 

V := + z//,2, 

and this definition does not depend on the choice of the complete energy orthonormal 
system {h\,h2\. By construction all energy measures Vj^g are absolutely continuous with 
respect to v and have integrable densities r(/, g). In particular, we can find Borel versions 
Zij of the functions r(/ij, /ij) and a Borel set C K such that for any x G Kq, the real 
(2 X 2)-matrix 

'Z'x '■= {Zij(x))ij =1^2 

is symmetric, nonnegative definite and of rank one, and we have ulK \ Kq) = 0. For 
X G K\Ko we may define to be the zero matrix. See for instance [IHl [651 [98] . Note that 
every Z^, x & K, acts as a projection in M^, and for fixed x the space (M^/fcerZ^., (■, Zrc-)^2) 
is isometrically isomorphic to the image space (^^(IR^), {■,Zx-)^2)- In addition, we may 
assume Kq is such that h{x) 7^ for all x G Kq. 

According to [521 Theorem 2.4.1] the regular resistance form {S,J^) defines a local 
regular Dirichlet form on L2{K,v). Therefore the measures Uf, f & J^, coincide with 
the energy measures in the sense of L33J, and the operation {f,g) 1— )■ T{f,g) taking two 
members f,g & into the density r(/, g) coincides with the carre du champ in the sense 

of [ig. 

Remark 2.1. Here we consider the L2-space L2{K, v) with respect to the Kusuoka measure. 
Note that the energy measures vj^g are singular with respect to the naturally associated 
renormalized self-similar Hausdorff measure on K, see [THl 1^ [15]. 

Setting 

h{x) := {hi{x) , h2{x)) and y := h{x) 

we obtain a homeomorphism h from K onto its image h{K) in M^, and the latter may be 
viewed with coordinates y. The collection of functions of form f = Foh with F G C^(M^) 
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is dense in J-", and for any such function the Kusuoka-Kigami formula 



(5) 



\ZVF{y)\l,dv 



K 



holds, where V-F is the usual gradient of F in M^. More generally, by the chain rule 
Theorem 3.2.2] the energy measure Vf of such / is given by \Z'W F {ii)\-^2dv . 

By L2,c(-^; ^) we denote the natural complexification L2(X, v) + iL2{X, u) of L2{X, v 
The closed form E on L^iX^ v) can be complexified by setting 



(6) 



■= £ifi,9i) -i£ifi,g2) + i£if2,gi) + £(91,92) 



for any f = fi + if2 and 9 = 9i + i92 from J^c '■= + ■ This yields a positive 
definite quadratic form £. on L2,c(-^; ^\ That is, £■ is conjugate symmetric, linear in the 
first argument, and £(j) > for any / G J-ic- We will use a similar terminology for the 
mappings considered in what follows. The form £■ is densely defined and closed. Similarly, 
and in a way consistent with ([6]), also the energy measure v^^g and their densities r{f,g) 
can be complexified. 

Consider J-c ® J^c endowed with the symmetric bilinear form 



(7) 



(a (g) 6, c (i) 



n 



bd r(a, cjdu, 



K 



a®b,c® d E J-c ® J^c- Let "He denote the Hilbert space obtained by factoring out trivial 
elements and completing. Following [191 120] we refer to it as the space of 1-forms on K. 
For simplicity we will not distinguish between an element a ® 6 and its equivalence class 
in "He- 

To rewrite several items in coordinates we also define the space 

Sc := span [f ® g : f = F o h, g = G o h with F, G G C^lM^)} . 
Theorem 2.1. The space Sc is dense in 'He- 
Proof. Note first that the collection Sq of elements f ® g with f = F o h, F E C,^(M^) 
and g G J-c is dense in l-Lc'- By the definition of He it suffices to approximate finite linear 
combinations ai®bi G J^c^J^c by elements of Sc- For fixed i let {fl"^)n be a sequence 
£^-converging to Oj. Then 

which converges to zero by the boundedness of the functions 6j. On the other hand every 
element '^j ® ^« of '^c can be approximated by elements of Sc'- For fixed i let {gi^^)n 
8- converge to 6j. The estimate ([3]) implies uniform convergence, and therefore also 

2 

^ (g) 6i - ^ (g) 51^"'' 

□ 



6,- 



5^ / 6.6,r(a, - /("V^^, 
"'-ft' 



/ {h - 9t''Y ^{ai,aj)du 
Jk 



goes to zero. 

Recall that we use the coordinate notation y 



y{x) = {hi{x),h2{x)). 
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Theorem 2.2. There are a family of Hilbert spaces {'Hc,x}x£x '^^^ surjective linear maps 
u ^ from He onto such that the direct integral V-cx^idx) is isometrically 
isomorphic to Tic and in particular, 

r 

11^11^= / \\^x\\n^'^{dx), uel-Lc- 
Jk 

For v-a.e. x & X the fiber Tic.x is isomorphic to C'^/ker Z^, and for any f = F o k and 
g = Goh with F,G e C^(M2) we have 

(8) \\f®9fn= [ \ZGiy)VF{y)\l,du. 

Jk 

Identity ([8]) obviously extends ([5]). 

Proof. For fixed x G Kq define a linear map (f^'-^c^ C'^/ker by 

Mg^f) ■■= Z,G{y)VF{y). 

Since h{x) ^ and linear functions belong to Cc(ffi^), the linear map (fx is surjective. 
Consequently there exists an isomorphism from Sc/ker ip^ onto C"^ /ker Z^, given by 

(9) ^x{U®9)x) = zMy)^ny). 

where (/ ® g)x = t^xif ® g) denotes the equivalence class mod ker (j)^ oi f ® g and the 
canonical epimorphism. From ^ we obtain 

(10) {f<^g)., = g{x){f®l)x 
for any f,gEJ^c and x E Kq. We write 

^c,x := S/ker (p^ 
and endow this space with the norm 

(11) \\if®gU^^^:=\ZxG{y)VFiy)y. 

Then becomes a isometric isomorphism between Hilbert spaces. For x E K \ Kq we 
set 1-Lic,x = {0} and k^. := 0. For every x E X the fiber Hc,x' is finite dimensional and 
therefore Hx extends uniquely to a surjective bounded linear map '■ "He 'Hc,x- For 
/ ® (? G iSc as above we have 



\\f®grH= / \gmf)du= / \G{y)\'\ZxVF{y)\iu{dx)= / \\{f ® g)x\\nxHdx). 
Jk Jk Jk 

Using bilinearity and the denseness of 5c, this extends to an isometric embedding 

k:Hc^ / UcxJ^idx) 
Jk 

of Tic into Hcx^idx). In fact k is onto: Assume that u E T-ic,xi^{dx) is such that 



= {ujj g)^= / {ujx,{f ® g)x)nx'^idx) = / g{x) {ujx,{f ®l)x)nx'^idx) 

J K ' J K 

for a\\ f ® g E Sc. Note that we have used ffTOl) . Then necessarily {ux, (/ Cg) '^)x)nx ~ ^ 
for Therefore ujx must be zero on 7^c,x for such x and integrating, we have w = 

in f^V-c,xi^{dx). □ 



8 MICHAEL HINZ^ AND ALEXANDER TEPLYAEV^ 

The definitions 

(12) c{a ^ b) := (ca) ® b - c{bd) and {a ®b)d := a ® {bd) 

for a,b,c G J-ic ® J^c and d G Loo,c(-^7 ^) extend continuously to uniformly bounded 
actions on Tic, 

(13) \M\n<ML^(x,m)\Mn and ||wc||^ < ||c||^^(^^„) , uenc- 
For w G "He and c G J-c we have the equality cu = uc in "He See [HI |20l [50]. By 

z/^,,,(A) := {ulA,r])H 

for uj,ri E "He and Borel set A C if, we define (weighted) energy measures for 1-forms. 
Note that for any / G J-ic we have i^/^i = z//. Note also that every u^^rj is absolutely 
continuous with respect to u, and x i— (w^;, ry^;)^ x ^ version of its density. To stress the 
similarity to the energy density we also use the notation 

(14) r^,x(Wa., T]x) ■= {(^x, Vx)'H,x ■ 

A derivation operator d : F<c — )■ T-ic can be defined by setting 

df:=f® 1. 

It satisfies the Leibniz rule, 

(15) dU9) = fdg + gdf, f,geJ^c. 
The linear operator d is bounded, more precisely, 

(16) lia/||^ = ^(/), /Gj-c. 

By the closedness of {S,J-'c) in L2^c{K,v) the derivation d may be viewed as an un- 
bounded closed operator from L2{K, v) into %. In coordinates y = y{x), the operator d 
agrees with the usual gradient operator V in M^. This can be phrased as a Corollary of 
Theorem 12.21 in terms of the isomorphisms from ([9]). It may be viewed as a 'pointwise 
formula' for the derivation d. 

Corollary 2.1. For any f = F oh and g = G oh with F G Cc(M^) and G hounded Borel 
measurable on we have 

^x[{gdf)x) = ZxG{y)VF{y) 

for u-a.e. x & K, and in particular, 

$,((9/),) = ZxVF{y). 

Because of the self-duality of "He we regard its elements also as vector fields and 5 as a 
generalization of the classical gradient operator. Let J-"^ denote the dual of J-c/ ~ with 
the norm 

=sup{|w(/)|:/G J-c, £{f)<l}. 
The space J^^ may be thought of as a 'space of distributions'. The symbol (■, ■) will denote 
the dual pairing between J^/ ~ and J-"^. Note that J-c C J^^ and {f^g) = {f,g)L2{Ku) 
for f,g e I'd ~. 
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Given a vector field of form gdf with f,g& Be, its divergence d*{gdf) can be defined 
similarly as in [47j by 

d*{gdf){^) :=- [ gT{f,if)du, 

J K 

seen as an element in J-"^. This extends continuously to a bounded linear operator d* 
from He into J^^ such that 

d*v{Lp) = - {v, d(f)^ , V e "He, f e Jb, 

see [m Lemma 3.1] for a proof. Note that here the divergence d*v is defined in a dzs- 
tributional sense. By restricting the space of test functions further this definition can be 
modified to fit into the context of distributions on p.c.f. fractals as studied in [83]. In 
coordinates we have the following expression. 

Theorem 2.3. Let f = F o h, g = G o h and u = U o h with F,U e C^(M2) and G 
bounded Borel measurable on M^. Then 

d*{gdf){u) = [ dwzAGVF)iU)iy)du, 

J K 

where for u-a.e. x & K , 

(17) div^,(GVF)(.)(y) := -$^%(a:)(GVF),(y)|^(y). 

( |T7|) may be seen as a bounded linear functional on C,^(R^). In a sense it remotely 
reminds of the divergence in Riemannian coordinates. 

Proof. We have 

-{gdf,du)^ = - [ {G{y)VF{y),Z,VU{y))^,u{dx) 

JK 

□ 

Let /S.y denote the energy Laplacian on K, that is, the infinitesimal generator of 
in L2{K,i'). Its complexification is the non-positive definite self-adjoint operator A,, 
uniquely associated with the closed form {£,J^c) on L2,c{K,i') by S{f,g) = — (/, A^^f) 
for f,gE J^c- Since —d*dg{f) = S{f,g) we observe 

(18) A,g = d*dg 

for any g G J-ic- Here f|T8l) is seen as an equality in J^^, below we will refer to an L2-context. 
In coordinates we have the following. 

Theorem 2.4. For any f = F ok with F G C^(R^) we have f G dam A^, and 

AJ{x) = TT{Z,}lessF){y) 
for u-a.e. x & K. Moreover, given arbitrary u = U ok with U G Cl^{E?), the identity 

f diYzA^F){U){y)du= [ AJudu 

J K J K 
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holds. 

Proof. The first statement is a special case of [981 Theorem 8]. For the second note that 
/ divz.(VF)(f7)(i/)z/(rfx) = - / {VF{y),Z,WU{y))^,uidx) 

J K J K 

A^fu dv. 

□ 

Remark 2.2. In the classical Euclidean case we have 

div(GVF) = GAF + VFVG 

for G G C^(M^) and F G C^(]R^), where V, A and div denote the Euchdean gradient, 
Laplacian and divergence, respectively. In our case we have 

d\gdf)=gA,f + VU\g) 

for f,g&J^, by [47l Lemma 3.2], seen as an equality in J^^. This may again be written in 
coordinates. Given f = F oh and g = G oh with F G G'^iM?) and G G G'^iM?) we have 

/ div^,(GVF)(f/)(y)rfz/ 
Jk 

= [ G{y)TT{Z,RessF){y)U{y)du + [ {VF{y), Z,VG{y))c2 U{y)du 
Jk Jk 

for any U G C^(]R^), as may be seen from the previous theorems. 

The divergence d* may also be seen in an L2-sense as an unbounded linear operator 
from into L2,c{K, v). As usual an element v G "He is said to be a member of dam d* 
if there is some v* G L2^(c{K,u) such that {v* , u) j^^^j^ ^-^ = — {v,du)y^ for all u G J^c- In 
this case we set d*v := v*. The operator —d* is the adjoint (codifferential) of d, and since 
9 is a closed operator, dam d* is dense in T-Lq- Note that for the previous distributional 
definition we obtain 

d*v{ip) = {d*v,ip)^^^j^^^y 

We end this section by a remark that allows to retrieve some more explicit information 
about the domain dam d* of the divergence operator d* in L2-sense. 

Remark 2.3. As the quadratic form {S, J^c) is closed, the operator d is closed. From (|3]) it 
follows that {£, J-'c) has a spectral gap. Therefore the range Imd of 9 is a closed subspace 
of "He and the space Tic decomposes orthogonally into Imd and its complement (Imd)'^. 
By orthogonality we observe (Imd)-^ = herd*. In |17] and |50] the space herd* had been 
identified as the space of harmonic forms (or vector fields) in the sense of Hodge theory 
if the topological dimension is one. 

Let dam Ai, denote the domain of A^, in L2,c{K, v). We obtain the following explicit 
description of dam d* . 
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Corollary 2.2. We have 

dom d* = {v E Tic : v = df + w : f E dom , w E ker d*} . 
For any v = df + w with f G dom and w G ker d* we have 

d*v = A,/, 

and for any g G dom Aj, formula [T^) holds in L2{K, u) . 

Note that the right hand side of this identity is explicitely seen to be dense in T-L: The 
denseness of the range of Green's operator in L2{K, p) can be used to see that dom A^ is 
dense in Jx- Therefore the image d{dom Aj^) of dom A,y under the derivation d is dense 
in Im d. 

Proof. According to Remark 12.31 any v G dom d* G H admits a unique orthogonal decom- 
position V = df + w with / G J-'c and w G ker d* . Since 

A J = d*df = d*v 

is in L2{K, v) we obtain / G dom A^. Conversely, any vector field v = df + w E l-L with 
/ G dom Ay and w G ker d* is a member of dom d* since for any u G J-'c we have 

-{v,du)^ = {df,du)^ = £{f,u) = -(A^/,m)^^(^^^). 

The last statements of the Corollary are obvious consequence. □ 

Remark 2.4. 

(i) Assume f = F o h. If F G C^(M^) then df G dom d* by Theorem [Ml For 
general F G C^(M^), however, we will not have / G domA^, and therefore also not 
df G dom d* . 

(ii) Even if / = F o k with F G C'^iM?) a simple vector field gdf can generally 
not be expected to be an element of dom d*. Let Vimd denote the orthogonal 
projection in Tic onto Im d. Given some measurable (or energy finite) g, the 
operator / h-j- Vimdigdf) is a complicated non-local first order pseudo-differential 
operator. See [l9j for some first results concerning pseudo-differential operators 
on fractals. 

3. DiRAC OPERATORS 

The definitions of the derivation d and the divergence d* allow to define related Dirac 
operators. Consider the Hilbert space 

and write 

{{f,uj),{g,v))^ := {f,g) 

for its scalar product. We define the Dirac operator associated with {£,J^) to be the 
unbounded linear operator D on L2,c{K, u) © "He with domain 

dom D : = J-'c © dom d* , 

given by 

D{f, Lo) := {-id*io, -idf), (/, u) G dom D. 
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In matrix notation we have 

Remark 3.1. Here the signs and imaginary factors are chosen in a way that fits with the 
following section. Of course other, possibly more physical choices, can be considered by 
obvious modifications. 

The next theorem was obtained in abstract form in [T9], but for the convenience of the 
reader we sketch a proof. Note that the results in Section [2] show that D is a pointwise 
defined local operator. 

Theorem 3.1. The operator {D, dam D) is self-adjoint operator on L2,c(-^? i^) © 'He- 
Proof. Obviously dam D is dense in L2^c{K, v) © "He- Recall that by definition 

dam D* = {{g, i]) E L2,c{X, m) © "He '■ there exists some ((?, 77)* G L2^<c{X, m) © "He 
such that ((/, w), {g, r])*)^ = {D{f, u), {g, r]))^ for all (/, u) G dam D} . 
Note that {—id)* = —id*. For arbitrary (/, and {g,ri) from dam D we have 
{{-td*uj, -idf), {g, T]))^ = {-td*u, g)L^(^x,m) + (-^^/' v)h 

= (w, -idg)^ + (/, -id*r])L,ix,m) 
= {{^J),{-id*ri,-idg))^. 
Consequently dam D C dam D* and D*u = Du for all u G dam D. □ 

Our next aim is to justify our nomenclature by showing that in an appropriate sense, 
:= D o D is the Laplacian. The Kusuoka Laplacian acts on functions, and we need 
to discuss a second Laplacian which acts on 1-forms. From [?8] we recall the following. 
Let 

dam 1 := {w G dam d* : d*UJ G J-c} 
and for oj G dom A,^ 1, define 

(19) A^,iw := dd*uj. 

To Aj, i we refer as the form Laplacian associated with {S,J^). The following result had 
been shown in jlHl Theorem 6.1]. For convenience, we briefiy sketch its proof. 

Theorem 3.2. Definition / flgj) yields a self-adjoint operator {A^^i, dom A^^i) on Tic- 
Proof. By Remark [2.31 we observe (Jm d)-^ C dom Aj, i. Let C be an £^-dense subspace of 
J-c such that for all / G C, we have A^/ G J-'c- Such a space always exists, for instance, 
we may use the image of J-'c under the Green operator A~^. For df with f E C and 
Ay/ = g E J-'c we have 

A,,,{df) = diAJ) = dg. 

Therefore also d{C) C dom A^^i. As C is £-dense in J^, its image d{C) is dense in Im d, 
and therefore dom Aj^ 1 is dense in "He- For any u G dom A^^i the identity 

(?7, Ay^iu;)^ = - ((9?7, (9*a;) ^^(^_^) = (A^^ir^, w)^ , 77 G rfom A^^i, 

showing dom A^ i C rfom A* and Al^u = A^^iu for all u G rfom A^ 1. □ 
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Now set dom := dom A^, © dom A^,^!, clearly a dense subspace of L2^c{K, v) © "He- 
Then the following is immediate. 

Lemma 3.1. For any (/, w) G dom ® dom we have D^{f,ijj) = {A^f, A^^ico) . 
In matrix notation this is 



D 



2 



A, 
A.i 



4. Magnetic Schrodinger operators 

The notions defined in Section |2] also allow to define a suitable generalization of the 
quantum mechanical Hamiltonian 

(20) H"^'^ = i-tV -Af + V 

for a particle moving in Euclidean space R" subject to a real [magnetic) vector potential 
A and a real valued (electric) potential V. The main result of this section is Theorem 14.11 
below, which tells that there exists a self-adjoint operator on L2^c{K, v) which generalizes 
fl20|) . Another result, Theorem 14.21 is a gauge invariance result and tells that, roughly 
speaking, we may restrict attention to divergence free vector potentials. 
We collect some preliminaries. Let 

'H = {v eUc-.v = v] 

be the space of real vector fields. An element a ^Ti may be seen as a bounded linear map- 
ping a : Fc — >■ "He defined by / i— )■ fa. I admits the estimate ||/ci||^ < i^) II'^II-h- 
For a ^Ti, define a bounded linear mapping a* : Tic — )■ J-J by f i— )■ a*v, where 



{a*v){ip) := I ipTnia^v) du = {aip,v)y^, e J'c- 
The boundedness follows from 



K 



sup 1 1 J ipT-^{a,v) du\ : (f e T^l ~ with E[}{>) — l| 
together with 



K 



^Yu{a,v) dv\ < \ML^(^K,u)\{(^^^)n \ <^(<^) 



1/2 



aWoj \\v 



where we have used ([3]) and the Cauchy-Schwarz inequality in "He. In terms of duality we 
observe a*v{ip) = {a*v,Tp). Given / G J-'c we have in particular a*af := a*{af) G -7-"^ with 



{a*af){(p) = / (pfTnia) du. 

J K 

Now let a G "H and V G Loo{K, u). The preceding considerations ensure that 
^"•''(/, g) ■■= ((-A. + '2^a*^ + a* a + i{d*a) + V)f, g) , 

f,g & J-'c, provides a well-defined form S""'^ on J-'c x J^c that is linear in the first and 
conjugate linear in the second argument. 

Proposition 4.1. Let a E H and V G Loo{K, u). 
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(i) For any f,g& J^c we have 

(21) ^'^■^(/, g) = {{-id - a)/, {-id - a)g)^ + {fV, g) . 

(ii) The form E""'^ is a quadratic form on J^c- 

Before proving Proposition I4.H we verify a product rule. For u G J^q and / G J^c define 
the product fu G by 

{fu){(p) := {u,ipf) , ipe Te- 
la particular, we ob obtain 

{f{d*a)){^) = {d*a,^) = -{a,d{^))^. 

Lemma 4.1. For a G "H and f G J-'c we have 

d*{fa) = f{d*a)+a*df, 

seen as an equality in J^^. 

Proof. Using the Leibniz rule (ITSl) we see that 

{f{d*a)){ip) = - (a,9(7V)>^ = - {a,lpdj)^ - (ajdip)^ = -{a*df){v) + d*{fa){v) 
for any Lp G J-'c- 

We prove Proposition 14. 1[ 
Proof. The right hand side of (12T|) rewrites 

(22) £{f,g) + M{f,g) + {fV,g), 
where 

M{f, g) := i {df, ag)^ - i {af, dg)^ + (/a, ga)^ . 
For the first summand on the right hand side we observe 



{idf, ag)^ = {df, -iag)^ = {-iag, df)^ 



and similarly for the second. Consequently M{g,f) = M{f,g), and and therefore the 
expression fl2^ is conjugate symmetric. We show that fl2^ equals £°''^{f,g). Note first 
that 

{af,dg)^= / fdTn{a,dg) = a*dg{f). 
Jx 

By Lemma [4.11 this equals 

d*{ga){f) -g{d*a{f) = -{ag,dj)^ - {d*a,gj)^ = - {df,ag)^ - {{d*a)f,g)^. 
Next, note that 

{df,ag)^ = [ gTn{df,a)du = a*{df){g) = {a*{df),g) . 

J K 

Consequently 

M(/, g) = 2i {a*df, g) + i {{d*a)f, g) + (a*a/, g) , 

and taking into account (JTSi) . identity (pTi) in (i) follows. Statement (ii) is a straightforward 
consequence. □ 
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Recall (dlD. To 

(23) U^-= [uj eH:v = v and r^,.(w) G L^{K, z/)} 

we refer as the space of (real) vector fields of hounded length. Assume a G "Hoo- Then the 
multiplication f ^ fa may be seen as bounded linear operator a : L2^c{K, v) — )■ He by 
/ I— J- /a, because 

(24) ||/a||^ = I \f\'TuAa)dv < \\^nM)\\L^^K,u) ■ 

From ( ITTi) it easily follows that 

5 = span{/(g)^: f = Foh,g = Goh with F, G G C^(R^)} . 
is a subset of Tioo- 

The estimate ( 12^ allows the application of classical perturbation arguments to prove 
the closedness of S""'^ and the essential self-adjointness of an analog of ( l20l) . 

Theorem 4.1. Let a G "Hoo o^^^^ ^ G L^oiK, v). 

(i) r/ie bilinear form (£^"'^, J-c) closed. 

(ii) r/ie unique self-adjoint non-negative definite operator associated with (£^"'^, J-^) 

i/"'^ = (-^5 - a)*(-i9 -a) + V, 
and the domain of A^, is a domain of essential self-adjointness for if"'^. 

Proof. Recall that £"''^{f,g) equals (122!) for any f,g E J-'c. We show that for any e > 
and any / G J-'c, 

(25) \M{f,f)\<e'S{f) + Ca,v\\f\\l^^^^), 

where Cay > is a constant bounded by plus a multiple of ||r^ .(a)||^^|,^^^. 

Then the result follows by the classical KLMN theorem, cf. [801 Theorem X.17]. By the 
boundedness of V clearly 

I {fV, f)L2(K,y) I — II^IIlooC/^,!') ll/llL2(i^,i/) • 

The bound covers the summand (/a, /a)^, and applying it once more, 

I (5/, I < \\dfU < WdfWn + ^ l|r«,-(«)llL.(x,.) II/IIL(;.,.))- 

Taking into account also f|T6l) . we arrive at fl25|) . □ 

Remark 4.1. Theorem 14.11 is a rather simple result and can certainly be improved. For 
instance, to obtain essential self-adjointness for magnetic Schrodinger operators on Eu- 
clidean spaces it is sufficient that the vector potential A is locally in L4 and div A is locally 
in L2, while the scalar potential V may be taken to be locally in L2. See [70]. What we 
would like to point out here is that any analog of such a hypothesis will again require the 
weighted energy measure Pa of a to be absolutely continuous with respect to the reference 
measure u. 
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If in fl20|) a gradient VA of a real-valued function A is added to the vector potential A, 
then the operator }j^+^^y is unitarily equivalent to H^'^ , more precisely, 

This property is usually referred to as gauge invariance, cf. [221 [70]. We observe a similar 
behaviour in our case. 

Theorem 4.2. Assume V G Loo{K, u) and a,b & 'Hoo- If b = a + d\ with some A G J-", 
then dom H^'^ = dom H"^'^ and 

^b,V ^ ^iX-^a,V^-iX_ 

Remark 4.2. Recall Remark 12.31 Let V± denote the orthogonal projection onto (Imd)-^. 
By Theorem 14. 21 the operator if"'^ is uniquely determined up to unitary equivalence by V 
and V±a. In this sense we can always restrict attention to divergence free vector potentials 
a G Hoo n ker d*. In the classical case the condition d*a = is referred to as Coulomb 
gauge condition, see for instance [87] . 



To prove Theorem 14.21 we first verify the following lemma. 
Lemma 4.2. Let a G Hoo, let X = Ao h, A ^ C^(M^) and set b := a + dX. Then we have 

£:<^'0(e-V) = ^''°(/), /e^c 

Proof. Note first that under the stated hypotheses we also have b G Hoo- 
Assume that f = Foh with F G C^{M.'^). By © we have 

$,,(-^9(e-*V)x) = -tZ,,i-tVAiy)e-'^(y^F{y)+e-'^^y^VF{y)) = $,((-e-^9A-ze-^"9/),) 

for i/-a.a. a; G X, and as the are isomorphisms also 

-^^{e'''f). = i-e-'^dX - le-'^df), 

in the spaces Hc^x- Integrating, we obtain 

-id{e"^f) = -e-'^dX - ie-'^df 

in He and therefore 

e'\-id - a)e-'^f = {-id - b)f, 
what implies the desired equality. For general / G J^c let (/n)n with /„ = F„ o and 
Fn G Cc(M^) be a sequence approximating / in S. Then 

lS-fi{f) - S^''{f^)\ = \S{f) - S{f^)\ + 2| {df, af)^ - a/„)^ | + | ||a/||^ - ||a/„||i |. 
The first term clearly tends to zero, for the second we have the estimate 

2\{df - dfn, af)^ I + 2| ((9/„, af - a/„)^ | 

< 28{f - f^fl' \\af\\^ + 2 11/ - fnW^^^K,.) Mn ^^P^ifY^^ 

11 

which tends to zero by ([3]). The third does not exceed 

\Hf - fn)\\nihf\\H+ 

what is similarly seen to converge to zero. □ 
We prove Theorem 14.21 
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Proof. For A as in the theorem we obviously have d\ G "Hoo- Let (a„) C iS be a sequence 
approximating a in "H and let {\n)n with A„ = A„ o /?, and A„ G C^(M^) be a sequence 
approximating A in Set 6„ := a„ + 9A„. Then Lemma [4.21 implies 

(26) g--y{e-'^"f)=S'-'''{f), fe7c. 
Now we first claim that 

(27) £:^'^(/)=lim^''-^(/), /GJ-c. 

n 

This follows from 

i^^v(j) _ ^^.,v^(j)| < 2| (5/, (h - 6„)/)« I + 1 ml - ii/^niii I 

because we have the upper bound 

2^(/)'/' ll/llL^(i.,.) \\b-hn\\n< 2^(/) II& - &nlk 
for the first summand, and for the second, 

ll(fe-MII (ll/&lk + ll/&n|| 

Combining, we arrive at (1271) . Next, note that we also have 

(28) £'''^(e-*V) = lim^'^"'^(e-^^"/), / G J'c- 

n 

To see this, note that 

|£"'^(e-^V) -^""'^(e"*^"/)| 

< |£"'^(e-^V) - ^'''^(e-^^"/)| + |£'''^(e-^^'7) - ^'^"•^(e"^^"/)! 

< 2| (9((e-^^ - e-^^")/), a>^ | + | ||ae-^V||^ - ||ae-^''7|| J, I 
+ 2| (9(e-^^"/), a - a„)^ | + | ||ae-^""/||^ - ||a.e-*""/||^ |, 

what can be estimated by 

(29) 2 \\d{l - e*(^-^"))e-'V||^ + 2 ||(1 - e'^^-^-'^)d{e~'^ f)\\^ 

+ 11(1 _ e'(^-^"))e-Va||^ + 
+ 2||9(e-^^"/)||^||a-a„lk 

+ ||e-^""/(a - a„)||^ (||e-*^'7a||^ + ||e"'''"/a||J • 
By the chain rule, [33l Theorem 3.2.2], we have 

^(l_e^(^-^")) = / |e'(^-^"f r(A-A„)rfz/ = ^(A- A„), 
Jk 

and consequently 

hm / |e^Vrr(l-e^(^-^"))dz/<lim||/||^^(^^,)£:(A-A„) = 0. 

The function 2; i— )■ 1 — e^*^ is continous and bounded, and by ([3]) we have lim„ A„ = A 
uniformly on K. Therefore 



lim / |l-e*(^-^")|2r(e-^^/)c/z/ = 
" Jk 
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by bounded convergence. Combining these estimates we see that the first hne in fl29|) 
converges to zero. The other terms in (129|) obey similar bounds, use ([3]), @ and note 
that 

sup£(e~*^") = sup / |e-'^"|2r(A„,)dz/ = sup£(A„) < oo. 

n n J X " 

They all converge to zero, and fl28|) becomes evident. Clipping (l26l) . (1271) and (l28l) we 
obtain the equality of closed forms 

which implies the coincidence of the associated self-adjoint operators and therefore The- 
orem 1121 □ 

5. Other fractals 

The results of the preceding sections easily carry over to regular resistance forms on 
finitely ramified fractals, see e.g. [SOI ISl ESI EB] for background and precise definitions. In 
this situation it is always possible to find a complete (up to constants) energy orthonormal 
system /ii, /i^ of harmonic functions. We assume that 

h{x) := {hi{x),...,hk{x)), x e X, 

defines a homeomorphism from X onto its image h{X) in M'^. We can then define a 
(normed) Kusuoka energy measure as the sum 

u := + ... + 

of the corresponding energy measures, cf. [98, Definition 3.5]. By [981 Theorem 3] the 
form {S,J^) is a regular Dirichlet form on L2(X, z/). Following Kusuoka [OS] one can 
construct a matrix valued function Z on X as before, such that Zij is a version of the 
density of i'hi,hj with respect to u. By [HSl Theorem 6] the collection of functions F o h 
with F G C^(M''') is dense in and formula ([5]) still holds. All results of Section 121 may 
be rewritten for X in place of the Sierpinski gasket K. The constructions of Sections 131 
and m do not depend on the specific structure of X and therefore remain valid as well. 

Remark 5.1. We would like to point out that by minor modifications Section HI applies to 
any local regular Dirichlet form that admits energy densities (i.e. a carre du champ) with 
respect to the given reference measure (see [151 HZ] ) • Section I3l does not require energy 
densities, although it is most naturally applicable for topologically one-dimensional spaces 
of arbitrary large spectral and Hausdorff dimensions (see [IS])- 
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